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Abstract. We introduce a method to relate a possible truncation of the star cluster mass function at the high mass 
end to the shape of the cluster luminosity function (LF). We compare the observed LFs of five galaxies containing 
young star clusters with synthetic cluster population models with varying initial conditions. The LF of the SMC, 
the LMC and NGC 5236 are characterized by a power-law behavior N d L oc L~ a dL, with a mean exponent of 
< a >= 2.0 ±0.2. This can be explained by a cluster population formed with a constant cluster formation rate, in 
which the maximum cluster mass per logarithmic age bin is determined by the size-of-sample effect and therefore 
increases with log(age/yr). The LFs of NGC 6946 and M51 are better described by a double power-law distribution 
or a Schechter function. When a cluster population has a mass function that is truncated below the limit given by 
the size-of-sample effect, the total LF shows a bend at the magnitude of the maximum mass, with the age of the 
oldest cluster in the population, typically a few Gyr due to disruption. For NGC 6946 and M51 this suggests a 
maximum mass of M max = 0.5 — lx 10 6 Mq, although the bend is only a 1-2 a detection. Faint-ward of the bend 
the LF has the same slope as the underlying initial cluster mass function and bright-ward of the bend it is steeper. 
This behavior can be well explained by our population model. We compare our results with the only other galaxy 
for which a bend in the LF has been observed, the “Antennae” galaxies (NGC 4038/4039). There the bend occurs 
brighter than in NGC 6946 and M51, corresponding to a maximum cluster mass of M max = 1.3 — 2.5 x 10 6 Mq. 
Hence, if the maximum cluster mass has a physical limit, then it can vary between different galaxies. The fact 
that we only observe this bend in the LF in the “Antennae” galaxies, NGC 6946 and M51 is because there are 
enough clusters available to reach the limit. In other galaxies there might be a physical limit as well, but the 
number of clusters formed or observed is so low, that the LF is not sampled up to the luminosity of the bend. 
The LF can then be approximated with a single power-law distribution, with an index similar to the initial mass 
function index. 
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1. Introduction 

The study of young extra-galactic star clusters has be¬ 
come a whole new field of research since the discovery 
of young massive clusters. The Hubble Space Telescope 
(HST) has made it possible to resolve these objects 
and to undertake systematic studies of the nature of 
these objects. Young clusters with masses in the range 
of our Milky Way globular clusters (10 4 — 1O 6 M 0 ) have 
been found in merging gala xies like the “Antennae” 

interacting galaxies like 
M51 Ternse i feOOOt iBastian et alJl2005bh . starburst galax- 
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ies iMeurer et al • US Ge GriLet ~ai]l2003ah _and even 
non-interacting spiral galaxies llLarsen fc R.ichtlerlll999f) . 

Recently, a relatively young and very massive clus- 
ter was discovered in the merger remnant NGC 7252 
dWhitmore et al.l Ill993l ). Its (dyn amical) mass was con - 
firmed to be as high as ~ 1O 8 M 0 llMaraston et alihoflll) . 
This suggests that there are clusters that fill the gap be¬ 
tween the mass range of star clusters and that of dwarf 
galaxies. However, it remains to be seen if every galaxy is 
able to produce such massive cluster or if there are phys¬ 
ical limitations to the maximum mass of star clusters in 
different environments. 
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The LF of star clusters is a powerful tool when study¬ 
ing populations of star clusters. It indirectly gives us 
information ab out the underlying mass function (MF). 
IZhang fc Fall ill 9991) showed for the young clusters in the 
“Antennae” galaxies that the cluster initial mass function 
(CIMF) can be well approximated by a power-law distri¬ 
bution: N(M) d M oc M~ a dM, with an exponent 1 of 
— a' = —2. They derived the ages of the clusters using 
reddening-free parameters, which is necessary to extract 
the initial mass function from the total luminosity func¬ 
tion. The exponent of the CIMF (— a') is found to be 
close to —2 in a wide r ange of galaxy enviro nments down 
to masses of ~ 10 3 Mq ( de Griis et all2003blh This resem - 
bles t he mass function of molecular clouds llSolomon et al.l 
ll 987th consistent with clusters forming from molecular 
clouds. The LF can also be approximated with a power- 
law distribution: N{L) d L oc L~ a d L, where values for 
the exponent —a are found i n a range of —2.4 up till —1.7 
ijT;arsenll20f)2l hereafter L02: lwiiituiord[2003lh In general 
the indices of the bright LFs are smaller (i.e. steeper) than 
the index of the CIMF (L02). One of the unanswered ques¬ 
tions is whether the difference in slope between the mass 
function and LF have a physical meaning or if it is a result 
of different measurement techniques and artifacts like the 
contamination by stars at lower luminosities. 

(Whitmore et al ] lll999fr observe a distinct bend in the 
LF of young star clusters in the “Antennae” galaxies, 
where the slope faint-ward is shallower than the slope 
bright-ward of the bend. The exact slopes depend on 
the different ways of correcting for stellar contamination. 
They argue that this could be the progenitor turn-over 
of the peak that is observed in the luminosity function of 
old globular clusters which appears at My ~ — 7.2, 

corresponding to a mass of 2 x 10 5 Mq. The difficulty 
with directly relating the LF to a CIMF is that the 
LF consists of clusters of different ages, and the mass- 
to-light ratio of clusters changes drastically when clus¬ 
ters age. Between 10' and 10 9 year, a star cluster of 
constant mass fades ab out 4 magnitudes in the F-band. 
Recently, iMengel et al.1 1 2005lf observed the clusters in the 
“Antennae” galaxies in the Ks-band. They also find a dou¬ 
ble power-law behavior in the LF and argue that the mass 
function has a turn-over. 

When there is no physical limit to the cluster mass or 
if there are not enough clusters to sample the CIMF up to 
any such limit, the mass of the most massive cluster will be 
determined by t he total n umber of clusters and the slope 
of the CIMF I Hunter et alJ 20031 her eafter H03). A simi¬ 
lar idea was posed by Weidner et al ] ll2004 . who suggest 
that the maximum cluster mass in a galaxy is determined 
by the star formation rate in that galaxy. It has not been 
shown yet in a large sample of galaxies that the most mas¬ 
sive cluster is a result of size-of-sample effects. H03 only 
showed that this is the case in the SMC and the LMC. In 


this study we add M51 to the sample in order to see if the 
most massive object in this galaxy is also a result of size- 
of-sample effect or if there is a physical limit above which 
clusters cannot form in this galaxy. In addition, we use a 
method to relate a possible truncation of the mass func¬ 
tion at the high mass end with the shape of cluster LF. 
To this end we introduce an analytical model to generate 
a cluster population and derive the LF for different clus¬ 
ter formation histories, disruption mechanism and CIMFs. 
We will show that even if the initial mass function of clus¬ 
ters is truncated, the brightest clusters in the sample may 
will sti ll be d eter mined by sampling statistics as found by 
Iwhitmorel 1 20031) and L02. 

The paper is organized as follows: In § 0 we intro¬ 
duce the data used for this study. In § El we present the 
masses and luminosities as a function of log(age/yr) for 
three galaxies in the sample. The LFs of cluster popula¬ 
tions in five different galaxies are presented in § El We 
introduce in §0 a cluster population model with which we 
can reproduce LFs based on various variables. The rela¬ 
tion between maximum cluster mass and maximum cluster 
luminosity is discussed in § El A discussion in presented 
in § [3 The conclusions are presented in § El 


2. Data 

2.1. SMC and LMC 

The ages and masses of clusters in the LMC and SMC were 
derived by H03 from ground based broad-band photome¬ 
try. The dat a was t aken by the Michigan Curtis Schmidt 
telescope bv lMassevl ( 2002 1 and covers a large field of the 
LMC and SMC (11.0 kpc 2 and 8.3 kpc 2 , respectively). 
Aperture photometry was used to derive U,B,V and R 
magnitudes. The colors and magnitudes were corrected 
for foreground reddening ( E(B — V) = 0.13 mag for the 
LMC and E(B — V) = 0.09 mag for the SMC). The ages 
an d masses w h ere derive d by comparing the colors with 
the lLeitherer et ahl ( 1999 1 models. For more details on the 
data reduction and age derivation we refer to H03. 

2.2. NGC 5236 and NGC 6946 

Data for these two galaxies was taken from L02 and we 
refer to that paper for details about the data reduction 
and analysis. Briefly, clusters were identified on F-band 
equivalent (F555W and F606W) HST/WFPC2 images 
and aperture photometry was carried out using a 3 pixels 
aperture. Cluster candidates were confirmed by a visual 
inspection of the images. For NGC 5236 and NGC 6946 
there were 1 and 3 HST/WFPC2 pointings available, cor¬ 
responding to 7.1 kpc 2 /53 kpc 2 , respectively. 

2.3. M51 


1 Through the remainder of this work we will use a’ for the 

exponent of the cluster initial mass function and a for the 

exponent of the total cluster luminosity function. 


We used the age, mass and luminosity d ata of m ore 
than 1 000 clusters in M51 found by iBastian et alJ 
(l2005hlf . Briefly summarized, aperture photometry was 
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performed on F336W, F439W, F555W, F675W, F814W, 
F110W and F160W (roughly [/, B,V, R, /, J and H) im¬ 
ages of the HST/WFPC2 and HST/NICMOS3 cam¬ 
eras. The magnitudes and colors were compared 
with the GALEV simple stellar population models 

1 Anders fc Fritze-v. Alvenslebenll2003t ISchulz et al.lf2002l) 

for different combinations of age and extinction. For de¬ 
tailed tests of the accuracy of the ages and the depen- 
de nce on photometric er rors and choice of models we refer 
to lBastian et al.l ll2005bh . With 2 HST/WFPC2 pointings 
we cover 69 kpc 2 of the inner region of the galaxy. 

3. Masses and luminosities at different ages 

3.1. Maximum masses at different ages 

3.1.1. Predictions 

The most direct way to study a possible truncation of 
the cluster mass function at the high mass end is to 
look at the masses of clusters as a function of the log¬ 
arithm of the ages. When the cluster initial mass function 
(CIMF) is a power-law distribution (NdM oc M~ a d M) 
and not truncated, the maximum mass in the sample is 
purely determined by the size-of-sample effect. A single 
age distribution of clusters with a minimum mass M m ; n 
consisting of N clusters, will statistically have a maxi¬ 
mum mass of M max = N 1 ^ 01 _1 ) M m ; n . H03 showed that 
the maximum mass per logarithmic age bin should in¬ 
crease with log(age/yr) (assuming a constant cluster for¬ 
mation rate), since each logarithmic age bin corresponds 
to a longer time interval at older ages and hence more 
clusters. The number of clusters at higher log(age/yr) in¬ 
creases when assuming a constant cluster formation rate 
as: ATdlog(f) oc tdlog(f). So, the maximum mass of a 
sample of clusters formed with a constant cluster forma¬ 
tion rate should increase with age (t) as 

Af max dlog(f) oc f 1/(a -1) dlog(f). (1) 

This is illustrated in Fig. El We statistically generated 
10 000 clusters from a CIMF with exponent —2 between 
6 < log(age/yr) < 10. In the left panel we plotted the 
results as a function of age and in the right-hand panel as 
a function of log(age/yr). The increase of the maximum 
mass predicted by Eq. ^is only visible in the right-hand 
panel. 

3.1.2. Observations 

For three galaxies in our sample (SMC, LMC and M51) we 
have ages and masses of the individual clusters available 
from earlier studies. For details on the derivation of the 
ages we refer to § Q and references therein. In Fig. O we 
plot the masses of the clusters as a function of log(age/yr) 
for the LMC (left), SMC (middle) and M51 (right). The 
completeness limit increases with age to higher masses 
due to fading of clusters caused by stellar evolution. The 
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Fig. 1. Result of a Monte Carlo simulation where 10 000 
clusters are sampled from a power-law mass distribution 
and evenly distributed between 0 < t < 10 10 . In the 
left-hand panel the masses are plotted vs. age and in 
the right-hand panel vs. log(age/yr). In the latter bi- 
logarithmic plot an increase is visible with a slope of 1 
(dashed line), which follows from Eq. |U 

first impression may be that the unequal number of clus¬ 
ters in each age bin is incompatible with a constant CFR. 
However, the observed age distribution is a function not 
only of the cluster formation rate, but also of the de¬ 
tection limit and disruption of clusters. In addition, the 
age-fitting method yields some we irregularities. For ex¬ 
amples all galaxies seem to have an over-density of clus¬ 
ters at l og(age/y r)=7.2. This is a known fitting artifact 
(e.g. Icieles et a,lJl2005a ). Therefore, the unequal number 
of clusters does not contradict a constant cluster forma¬ 
tion rate. 

The dashed line shows the slope of the predicted in¬ 
crease of the maximum mass, based on a' = 2 and Eq. [TJ 
For the LMC and SMC the maximum cluster mass as a 
function of log(age/yr) can be well approximated by the 
size-of-sample prediction, as was shown already by H03. 
For M51 however, we see that for ages higher than 10 8 
year, there are more old high mass clusters predicted based 
on the size-of-sample effect than observed. However, the 
observations show that there are no clusters more massive 
than ~ 10 6 Mg. 

We also performed a linear fit to the maximum mass 
vs. log(age/yr). This is shown in Fig. [3] The most massive 
cluster in every log(age/yr) bin is indicated with a dot. 
The fit is shown with a full line. The LMC clusters older 
than log(age/yr)= 9.5 were not included in the fit, because 
there seems to be an age gap in the cluster distribution be¬ 
tween 3 and 10 Gyr. For the SMC and LMC the fit is con¬ 
sistent with the size of sample effect and a CIMF exponent 
of — a 1 = — 2.4, which was found already by H03 for the 
SMC and LMC. This is a rather steep slope for the CIMF 
and H03 showed that a fit to the first bins (log(age/yr) 
< 8) yields a value of —a! = —2. It might be that for 
higher ages we are seeing a lack of old massive clusters. 
Note also that we here assume a constant formation rate 

















4 


M. Gieles et al.: The luminosity function of young star clusters 





Fig. 2. Cluster masses for different log(age/yr) for the LMC (left), the SMC (middle) and M51 (right). All samples 
are limited by a detection limit. Since clusters fade when they age, the lower mass limit goes up in time, the expected 
increase of the maximum cluster mass as a function of time based on the size-of-sample effect and an exponents for 
the initial cluster mass function of a' = — 2 . 


of cluster and no mass loss of the clusters after formation. 
We will look at the effect of this assumption in § ra For 
M51 we find a much shallower slope (0.26) than expected 
on account of the size-of-sample effect. If this slope was 
caused by the mass function, it should have an exponent 
of —a 1 = —4.3, which would be a much steeper slope than 
found in any other study . In fact, the CIMF has been de¬ 
termined for M51 bv lBik et alJ 1 2003 1 from young clusters 
and was found to be —2.1 ± 0.3. 


3.1.3. The effect of mass loss, infant mortality rate 
and variable cluster formation rates 

In the conversion of the M max vs. log(age/yr) relation to 
a CIMF slope we assume that the clusters have not lost 
any mass after their formation and that the clusters were 
formed with a constant formation rate. In reality, only 10- 
30% of the clusters will survive the initial 10 Myr, clusters 
will lose mass due to stellar evolution (SEV) and tidal ef¬ 
fects and the cluster formation rate (CFR) might not be 
constant at all. To quantify these effect, we generated a se¬ 
ries of maximum masses consistent with the size-of-sample 
effect and a mass function exponent of — 2 , assuming a con¬ 
stant CFR. We then evolve each cluster mass as a function 
of i ts age accordi n g to th e analytical formulas introduced 
by ILamers et alJ ( 2005a .'). The mass of the cluster as a 
function of time can be well approximated by 


M p (t) 


(M t/4ev F- 7l0 47 



( 2 ) 


where M p (t) is the present mass of the cluster as a func¬ 
tion of its age, M ?: is the initial mass of the cluster, 
/x sev = M(t)/Mi is the fraction of remaining mass after 
mass loss due to stellar evolution has been subtracted. The 
value of fj, sev is based o n the GALEV models and analyti¬ 
cally parameterized (see lLamers et al ] ( 2005al) for details). 


ti is the disruption time of a 10 4 Mq cluster and 7 is a di¬ 
mensionless index indicating the dependence of the disrup¬ 
tion time of clusters on the initial mass; tdis oc M?. The 
value for 7 is found to be_0.62 based on_observations and 
Wbody simulations I Lamers et al.l2005bl) . EquationElde- 
scribes the mass of a cluster as a function of time based 
on a variable final disruption time, taking mass loss due 
to stellar evolution and tidal evaporation into account. 

We consider two scenarios: 1.) the SMC/LMC case and 
2.) the M51 case. This is because the masses and dis¬ 
ruption times are very different. For the SMC/LMC case 
we choose the maximum mass to be M max = 10 3 Mq at 
log(age/yr) = 6.5, corresponding to the observed maxi¬ 
mum mass at young ages. For the M51 case we choose a 
higher maximum mass; M max = 5 x 10 4 Mq at log(age/yr) 
= 6.5, also based on the observed masses at young ages. 
We then assume a typical disruption time of clusters for 
the SMC and LMC of £4 = 3 x 10 9 yr and Lj = 2 x 10 s 
yr for M51. The disruption times are based on the o b- 
servat ionally dete rmined value s of iBo utloukos fc Lamer! 
ll2003l) (SMC) andlGjeles_eCalJ |2005a|) ( M51) and the pre- 
dicted values of iBaumgardt fc Makinol ll2003l) (LMC). In 
Fig. El we show the slope of the cluster maximum masses 
as a function of log(age/yr) as predicted by the size-of- 
sample effect with a full line and the shallower increase 
after the clusters have been evolved in time for different 
scenarios. In the top panels we show the effect of stellar 
evolution and disruption for the two different scenarios. 
The disruption time of clusters in M51 is almost a fac¬ 
tor of 15 shorter than for clusters in the SMC/LMC, but 
since the masses are higher and disruption is mass depen¬ 
dent (fdis ex M 9 - 62 ), the effect of the evolution is similar in 
both cases. Since this simple model can easily be extended 
to more realistic scenarios, we now include the effect of a 
high infa nt mo r tality ra t e in M51 in the bottom left panel 
of Fig. El iBast.ian et al . (l2005hl) estimate that 70% of the 
clusters does not survive the first 10 Myr, independent 
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log(binsize/Myr) 

-0.2 0.8 1.8 2.8 3.8 



Fig. 3. The most massive cluster (dots) in log(age/yr) bins 
of 0.25 clex for all three galaxies. The fit to the most mas¬ 
sive clusters is shown with the full line in all panels. The 
values of the slopes are indicated. The top panel shows 
how the binsize increases with increasing log(age/yr), in¬ 
troducing the size-of-sample effect. 


of their mass, so we increase the maximum mass for bins 
below 10 Myr with a factor of 3.3 to simulate the higher 
masses at young ages. The result is shown in the bottom 
left panel of Fig.0 The inclusion of infant mortality makes 
the slope again slightly shallower. However, still not even 
close to the observed value. The last effect we consider, is 
a steady increase in the CFR. Since the maximum mass 
is proportional to the CFR, an increase in the CFR could 
result in more young massive clusters. To get to the slope 
of 0.26 as observed in M51 (bottom panel of Fig. 0 , we 
apply an additional increase in the CFR. We assume that 
increase exponentially in time as CFR oc t s . To get to a 
slope of 0.26 as observed, we need 5 = —0.5, which means 
the CFR has increased a factor of 22 in the last few Gyr. 
M51 is in interaction with NGC 5195 which might have 
enhanced the CFR. In earlier work, however, we have de¬ 
termined the CFR increase to be in the order of a factor of 


3-5 ijGieles et alll2005ali. which is in good agreement with 
independent determinations of the increase of the star for¬ 
mation rate for these kind of interacting s ystems, based 
on Ha measurements feergvall et al • m3). Also, the in¬ 
crease in the CFR seems to be more in bursts rather than 
an exponential increase in time. A burst in the CFR will 
not affect the slope of Fig.0 it will just show up as a peak. 
Worth noting is that the two peaks of M51 (Fig. 0 at 
log(age/yr)=7.8 and log(age/yr)=8.8. They coincide with 
the predicted moments o f encounter with NGC 519 5, ~75 
Myr and ~450 Myr ago ISalo & LaurikaineifeoQ Oj). 

With this we show that part of the deviation from a 
slope of 1 in Fig. 0 is caused by stellar evolution, cluster 
mass loss due to tidal effects and infant mortality rate, 
but even the short disruption time and high infant mor¬ 
tality rate of M51 can not explain the observed shallow 
increase of the maximum cluster mass with log(age/yr) 
(Fig.0 bottom), without assuming an unrealistically high 
increase in the CFR. 
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Fig. 4. A simulated set of initial maximum masses as a 
function of log(age/yr) (full line) for the SMC/LMC case 
and the M51 case with different effects. The full line rep¬ 
resent what is expected based on the size of sample effect 
and —a' = —2 and therefore has a slope of 1 (see Eq. 0. 
The initial masses were evolved in time taking into ac¬ 
count mass loss due to stellar evolution and tidal effects. 
The evolved masses (dots) are fit and the slope is indicated 
with dashed lines. 


3.2. The cluster mass function 

If the CIMF is truncated, the total MF of clusters should 
also show a truncation. As a check, we make the mass func¬ 
tion, i.e. dN/dM as a function of mass, for the SMC, LMC 
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Fig. 5. Mass functions of the LMC, SMC and M51 for an 
age and mass limited sample. The result of a truncated 
power-law fit is over-plotted and the corresponding M* 
is indicated. The xi of the truncated and untruncated 
power-law fits are compared by showing the ratio. 


and M51. Since the detection limit cuts the sample at dif¬ 
ferent masses (Fig. [2J, we need to generate a mass and age 
limited sample. We choose minimum mass and maximum 
age combinations of [log(M min /M©), log(age max /yr] = 
[3.1,9], [3.5, 9] and [4.2, 8.7] for the LMC, the SMC and 
M51 respectively. The resulting mass functions are shown 
in Fig. 0 The three mass functions where fitted with a 
single power-law, with variable slope and vertical scaling 
and with a power-law with exponential cut-off: NdM oc 
M~ a exp (—M/M*) d M, where the cut-off mas s is an 
extra variable. This function is similar to the ISchechteil 
dl97fih function, which will later be used to fit the lumi¬ 
nosity function of clusters. The xt results of the trun¬ 
cated and untruncated power-law fits are compared by 
showing the ratio. The fit of the truncated power-law is 
over-plotted. For the SMC and LMC, however, this fit is 
similar to an untruncated power-law fit. The M51 mass 
function is better fitted with a truncated power-law. The 
values we find for —a' are —2.05 ± 0.01, —2.03 ± 0.02 and 
— 1.70 ± 0.08 for the SMC, LMC and M51 respectively. 
The value for M51 can be explained by mass dependent 
cluster disruption. If star clusters are formed with a CIMF 
with index —2, then the total mass function will have the 
same index, as long as the disruption time scale is longer 
than the maximum age in the sample. For the SMC and 
LMC the maximum age in the sa mple is 10 9 yr, and the 
disruption time scale is >> 1 0 9 yr dHaumga.rdt Rr. Ma.kinol 
l2003tlBoutloukos fc Lamersll2003l) . For M51, however, the 
disruption time is shorter t han the maximum age in the 
sample dCieles et al.ll2005ai) . Since the low mass clusters 
disrupt faster than the high mass clusters, the mass func¬ 
tion gets shallower at older ages. Therefore we find that 
the MF has an index which is g reater than —2 .1, which is 
the index of the CIMF of M51 ( Bik et aHl2003 '). 


3.3. Maximum luminosities at different ages 

Another way of looking at the size-of-sample effect, is by 
plotting the My magnitudes of all clusters as a function 
of log(age/yr). In Fig. |3we plot the magnitudes for the 
three cluster populations shown in Fig. [3 The advan¬ 
tage of looking at magnitudes instead of masses, although 
they are directly linked by the age and mass dependent 
magnitudes from SSP models, is that this figure helps to 
understand the LF, which will be discussed in § □ The 
solid line indicates the fading of a single mass cluster 
based on the GALEV models, i.e. without tidal disrup¬ 
tion. We shifted the starting magnitude to the brightest 
clusters at log(age/yr) = 6.5. The dashed lines now repre¬ 
sent the luminosity of the most massive cluster as a func¬ 
tion of log(age/yr) based on the size-of-sample effect and 
a' = 2.0. From the GALEV simple stellar population we 
took the evolution of the magnitudes in time for a single 
mass cluster with some reference mass (M re f), which is 
the maximum mass in the age range log(age/yr) < 7. We 
then add extra mass to M re f as a function of log(age/yr), 
following Eq. ^ The magnitude then depends on the age 
and the exponent of the CIMF as 

M v (t) oc M v {M ieU t ) - 2.5 logtf 1 /^'- 1 )) (3) 


where M le f is a reference mass based on the observed max¬ 
imum mass at young ages, My(t ) is the maximum mag¬ 
nitude as a function of age, t is the cluster age and a' 
is the exponent of the CIMF. The dashed lines represent 
the maximum luminosity based on the size-of-sample ef¬ 
fect and Eq. |3 for a' = 2. The lower luminosity increases 
with log(age/yr) for M51. This is since the clusters of M51 
had to be detected in the F439W band in order to make 
the observed sample. Since clusters fade more rapidly in 
the filters bluer than the V band, we miss some clusters 
at older ages and faint magnitudes. Since we are here only 
interested in the maximum luminosity, this will not affect 
our results. When the colors are used for cluster identifica¬ 
tion, these color differences have to be taken in to account. 

The maximum luminosities of clusters in the LMC and 
SMC follow the dashed lines better than the fading line 
of the single mass cluster. The most luminous clusters in 
M51 seem to follow the fading of the single mass M re f 
better then the dashed lines which are based on the size- 
of-sample prediction, which suggests that the underlying 
maximum mass is truncated at that mass. 

Interestingly. IZhang fo FalJ ( 19991) made a similar plot 
for the luminosities of clusters in the “Antennae” galaxies 
(their Fig. 2). They also over-plot the luminosity evolution 
of single mass clusters in time. Their maximum luminosi¬ 
ties seem to follow the models very well, which suggests 
that the maximum luminosity and hence the maximum 
mass of clusters in the “Antennae” galaxies does not follow 
the predicted size-of-sample relation of H03. All this sug¬ 
gests that the maximum mass of clusters in M51 and the 
“Antennae”galaxies are not determined by size-of-sample 
effects, but instead has a physical upper limit. The most 
luminous cluster in these galaxies however, in general be- 
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Fig. 6. Cluster absolute magnitudes ( My ) for different ages for the LMC (left), the SMC (middle) and M51 (right). 
All samples are limited by a detection limit, which is why the lower parts of the diagrams are empty. The solid 
line shows the fading of a cluster in the F-band based on the GALEV models. The dashed line shows the expected 
luminosity of the maximum cluster mass as a function of time based on the size-of-sample effect and an exponents of 
the initial cluster mass function a' = 2. 


ing a young clusters due to the rapid fading in time of 
star clusters, is determined by size-of-sam pl e eff ects. This 
agrees with what was found bv lWhitmora (2 0031) and L02, 
who showed that the maximum cluster luminosity in a 
galaxy depends on the total number of clusters in the 
galaxy above a certain luminosity. We will continue with 
this subject in §0 


which are 1.5 mag brighter (My = —5 for the LMC and 
My = —6 for the SMC). H03 only corrected for foreground 
reddening ( E(B — V) = 0.09 and 0.13 for the SMC and 
LMC respectively) and we have no extinction estimates 
for the individual clusters available. The detection limit 
is shifted the same amount as the LF when correcting for 
foreground extinction. 


4. The luminosity functions of five galaxies 4.1.2. NGC 5236 


In this section we will present the LFs of (young) clus¬ 
ters in the SMC, LMC, NGC 5236, NGC 6946 and 
M51. For each galaxy the number of clusters as a func¬ 
tion of My was determined and a fit of the form 
N dMy = A 10^ Mv dMy was performed. The index /3 of 
the magnitude distribution relates to the exponent of the 
LF as 

-a = -2.5/3- 1 (4) 

where a is the exponent of the LF (N dL cx L~ a dL). In 
the following subsection we will first present the LFs and 
fit a power-law function to all cluster samples. In § m we 
will fit alternative functions and compare them with the 
power-law fits to determine the best fit model to the LFs. 
All results are summarized in Table 1. 

4.1. The observed luminosity functions 

4.1.1. SMC and LMC 

The LFs of the SMC and LMC are based on the data of 
H03. The LFs are presented in the top two panels of Fig. 0 
H03 report lower limit magnitudes of My = —3.5 and 
—4.5 for the LMC and SMC, respectively. Since the sample 
starts to be incomplete a bit bright-ward of the cut-off 
magnitude, we used very conservative completeness limits 


The cluster sample of NGC 5236 is presented in L02. From 
the completeness tests presented in that paper and a dis¬ 
tance modulus of 27.9, we found a 90% completeness limit 
of My = —5.5. We corrected the measured magnitudes 
for ga lactic foreground extinction based on lSchlegel et alJ 
( 1998 1. We used the sample which had spurious sources 
removed after visual inspection. A foreg round extinctio n 
correction of 0.218 mag was applied ( Schlegel et al . ml. 
The sharp cut-off at high luminosities is because some of 
the bright sources were saturated. 


4.1.3. NGC 6946 

The LF of clusters in NGC 6946 was presented in L02. 
We here use the sample of clusters where spurious sources 
have been removed after visual inspection. The 90% com¬ 
pleteness limit, following L02, is V = 22.0 and a distance 
modulus of 28.9 mag was used. The magnitudes were 
corrected for 1 .133 m ag foreground extinction based on 
ISchlegel et alJ (Il998l) . The final LF can be seen in Fig. 0 
(second from bottom). The LF contains more bright clus¬ 
ters than the previous ones. In addition the slope is slightly 
steeper. The fact that LFs with brighter clusters are better 
fitted with steeper power-laws was already noted by L02. 
In § 14.21 we will compare this fit with a double power-law 
fit and a Schechter function. 
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4.1.4. M51 

For the LF of M51 the data from iBastian et all ll2005bl) 
is used. We refer to this paper for detailed explanation 
on data reduction and source detection. Ideally, we select 
clusters for the LF independent of the applied age fit¬ 
ting routine. The best criterion for this would be to select 
on extended objects brighter than the 90% completeness 
limit in F555W. Since we only found reliable radius es¬ 
timates for 305 objects, we selecte d the observed cluster s 
that passed our fitting criteria from lBastian et all C 2005bh . 
e.g. the sources that were well fit by a cluster model. Since 
one of our source selection criteria is that the object has to 
be above the 90% completeness limit in at least four bands, 
it is of importance which band we use. Comparing the 
SEDs for different ages of the GALEV models to the 90% 
completeness limits in all filters, we found that all clusters 
are always above our 90% completeness limit in F555W, 
F675W and F814W once they are detected in the F439W 
filter. If we apply the F555W 90% cut to the sample, our 
four-filter criterion will cut the sample sometimes because 
the source has to be visible in F439W as well, depend¬ 
ing on the age. To be sure that for all ages our sample is 
limited by a cut in the F555W band only, we add the max¬ 
imum color difference F439W - F555W from the GALEV 
models for clusters between 4 x 10 6 and 10 9 years to our 
F555W 90% completeness limit. This makes the cut in this 
filter 0.6 mag brighter than the 90% completeness limit 
of 22.7 m ag. A distance modul us of M51 for 29.62 mag 
was used jFeldmeier et alj fl99?l) and a foreground extinc¬ 
tion correction of 0.117 mag was applied ijSchlegel et alJ 
1 998 :. The LF for all clusters brighter than 22.0 mag in 
F555W is shown in the bottom panel of Fig. 0 As was 
the case for the LF of NGC 6946, we also fit a power-law 
whic h is steep er than the power-law CIMF index —2 of 
Izhang fc Fall! 11 999 1. We will consider different functions 
in 8 14.21 

4.2. Fitting other functions to the luminosity functions 

The brighter LFs of NGC 6946 and M51 are fitted with 
power-law functions with steeper slopes (—a ~ —2.4) than 
the fainter ones (—a ~ —2). When a single power-law 
function is fitted to faint sides of the LFs of NGC 6946 
and M51, we find similar values as for the SMC, LMC and 
NGC 5236. This suggests that the LFs turn over at the 
bright side, something what was found for the “Antennae” 
galaxies by W99 as well. In order to quantify the deviation 
from a single power-law, we fitted all LFs of 8 14.1 1 with a 
double power-law function with variable break point and 
a Schechter function. The Schechter function is generally 
used for the LF of galaxies. It has a power-law nature 
on the faint end and a exponential drop at the bright 
side. In terms of magnitudes the Schechter function can 
be expressed as 


N(M\)dM\ = 0.4 ln(10) W10°- 4( "“ +1)(M *“ Mx) x 



M v mag 


Fig. 7. The (foreground extinction corrected) LFs for the 
five galaxies in our sample. In all graphs the dashed line 
represent the 90% completeness limit. The solid lines are 
fits to the data above the 90% completeness limit over the 
full range of the form NdMy = A If) 13 Mv dMy■ From this 
the slope of the LF can be found using Eq. Q] 


exp 


- 10 ' 


0.4 (M» —Mx) 


d M. 




(5) 


where M\ is the magnitude in a band with central wave¬ 
length A, AT* is a normalization number, a is the exponent 
of the power-law part of the LF and M* is the characteris¬ 
tic point of the function. M* is the equivalent of L* when 
using the Schechter function to fit the LF of galaxies. 

To compare the different function fits for all five 
galaxies, we compare the reduced y 2 values of the 
different fits. In the top panel of Fig. [ 8 ] we show 

the ratio of X^ s i n gl e p ower -law/X^ double power-law for the 
five LFs. In the bottom panel we show the ratio of 

single power—law/ "X-v, Schechter- The xl comparison shows 

that the LFs of NGC 6946 and M51 are better fitted with 
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a function that drops at the bright end, although from 
the values in Table 1 and the xt comparison mentioned, 
a single power-law distribution can not be excluded with 
a high significance level. 

In Fig. |21 we show the best fits to the LFs for all five 
galaxies, based on their xt results. The number of clusters 
is normalized to the observed area for fair comparison. We 
extrapolated the faint side of the LFs of NGC 6946 and 
M51 with a dashed line, this to show that these galaxies on 
the faint side have similar slopes as the other three cluster 
LFs. Another interesting fact is that NGC 6946 and M51 
are on top, due to a high CFR. Therefore, they are also 
the only galaxies with cluster brighter than My = —10 
mag. 




Fig. 9. The best fit result to all five LFs. The number of 
clusters is normalized to the observed area of the galax¬ 
ies in kpc 2 . The full line is the region where data points 
were available from observations. The dashed lines are ex¬ 
trapolations. The two galaxies which are better fitted by 
a double power-law (NGC 6946 and M51), are on top and 
the only galaxies with cluster brighter than My = —10 
mag. 



Mv mag 


Fig. 11. E( B — V") vs. luminosity for all clusters of M51, 
taken from lBastian et aD 1 20051:1 ). The average values at 
different magnitudes are indicated. 


4.2.1. A closer look at NGC 6946 and M51 
NGC 6946 


Fig. 8. Comparison between the xt results of the dif¬ 
ferent functions that were fitted to all LFs. Top: Ratio 

°f xt, single power-law 0Ver Xt, double power-law f° r the five 

galaxies of which we have LFs. Bottom: Ratio of 

xt, single power-law 0ver xl, Schechter' The LFs ° f NGC 6946 

and M51 are better fitted by a double power-law distri¬ 
bution or a Schechter function, as compared to a single 
power-law distribution. 


To get a realistic idea of the location of the breakpoint, 
which we can compare with other galaxies, we need to 
take into account possible reddening. L02 did not deter¬ 
mine reddening for the individual objects, but only men¬ 
tions the Galactic foreground reddening. For the direction 
of NGC 6946 this is Ay = 1.13 and we corrected the V 
magnitudes of each cluster for this amount. The break¬ 
point was determined by fitting double power-laws with 
a variable breakpoint, exponents and normalization con¬ 
stant. We found the break to be at My = —8.9 ± 0.4, with 
an exponent faint-ward of the break of —a± = —1.7 ± 0.5 
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Fig. 10. The extinction corrected LFs for NGC 6964 (top) and M51 (bottom). Results for a double power-law fit 
(left) and a Schechter function (right) to the same data are shown. 


and —02 = —2.4 ±0.1 bright-ward of the break. The (fore¬ 
ground) extinction corrected LF with the double power- 
law fit can be seen in the top left panel of Fig. m W99 
observe similar behavior for the LF of the “Antennae” 
galaxies. When they fit two power-laws, they find the 
breakpoint around My = —10.3, depending on the dif¬ 
ferent samples they selected._ 

We also fitted a Schechter function ( Schecht^ j1976 1. 
which has a more gradual bend and it is usually used to 
fit the LF of galaxies. This was also done for the clusters 
in the “Antennae” galaxies by W99. 

The top right panel of Fig. shows the result of 
the Schechter function fit to the (foreground) extinction 
corrected LF of NGC 6946. The bend in the Schechter 
function occurs at M* = —10.2 ±0.6, which is 1.2 mag 
fainter than the M* of the “Antennae” galaxies (M* = 
— 11.4 ±0.4). The break point in the double power-law 
distribution occurs 1.2 mag fainter as compared to the 
“Antennae”. 


M51 

For M51, we corrected all magnitudes for the extinction 
values w e acquired from the age fitting fsee lBastian et, all 
•: 200-51 for details). This is important, sinc e extinction ef - 
fects might also affect the shape of the LFs. lBastian et, ahl 


(l2005hlf found that young clusters have on average higher 
extinction values than the older ones. L02 looked at the 
effect of differential extinction on the LF. He found that 
it makes the LF steeper, but only when there is already 
a cut-off in the mass distribution on the high mass end. 
If the mass function is not truncated, differential extinc¬ 
tion does not affect the shape of the LF, since then still 
power-laws with the same index would be added, result¬ 
ing in a power-law with that index. In case there is a 
relation between the luminosity and extinction, for exam¬ 
ple because we can see more intrinsically bright clusters, 
the bright part of the LF will be extinct ed more and 
this could result in a drop. To check this, we pl otted the 
availab le extinction and luminosities from lBa,stia,n et alJ 
(' 2005bh in Fig. We averaged the extinction for differ¬ 
ent luminosity bins, and there is no relation between the 
two parameters. So the fits to the LF of M51 will have a 
higher reliability than the ones of NGC 6946, since there 
individual source extinction is not known. 

When we fit a double power-law distribution, the lo¬ 
cation of the bend is at My = —9.3 ± 0.4. The expo¬ 
nents on the faint and bright side of the breakpoint are 
—aq = —2.0 ± 0.1 and — 012 = —2.7 ± 0.3 respectively 
(bottom left panel of Fig. llOfl . These values are very simi¬ 
lar to those found for NGC 6946. The bottom right panel 
of Fig. shows the result of a Schechter function fit to 
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the M51 clusters. The exponent is —a = —1.9 and the 
bend is at M* = —10.3 ± 0.5. 

The results of all fits to the LFs are summarized in 
Table 1. We added the values found by W99 for the 
“Antennae” galaxies as a comparison. The two exponents 
of the double power-law fits are very similar between 
NGC 6946, M51 and the “Antennae”. The breakpoint in 
the case of the double power-law fit, however, occurs about 
1.2 mag brighter in the “Antennae” galaxies. The value for 
M* is also about 1 mag brighter compared to NGC 6946 
and M51. 

In the next section we will attempt to explain the 
shape of the LFs, in particular the bend observed for 
NGC 6946, M51 and the “Antennae” galaxies. 

5. A cluster population model 

Our goal is to understand why we observe a bend in the lu¬ 
minosity function of NGC 6946, M51 and the “Antennae” 
galaxies and not in the other galaxies in the sample. The 
luminosity function (LF) of a star cluster population is 
a combination of multi-age cluster initial mass functions, 
which all have age dependent mass to light ratios. Clusters 
will lose mass in time due to stellar evolution and tidal ef¬ 
fects, which will also affect the luminosity as a function of 
age. When there are no ages and masses of the individual 
clusters available, an alternative way to look at how the 
LF is built up is to model the LF based on various phys¬ 
ical input parameters. In this section we will analytically 
generate the LF of a star cluster population. 

5.1. Creating a synthetic cluster population 

5.1.1. Set-up 

The synthetic clust er populati on is based on an analytical 
model presented in iGieles et ali i'2003a i. In that work it 
was used to predict ages and masses of the cluster popula¬ 
tion of M51. The LF follows directly from these predictions 
since the combination of age and mass can be converted 
to luminosity using simple stellar population models. 

The model creates a series of cells equally spaced in 
log(f) and log(M). Then a CIMF and formation rate are 
chosen. Based on these choices each cell is assigned a 
weight (w(f, Mi)). To acquire a CIMF with exponent —a' 
and a constant formation rate, w(t, Mi) can be written as 

w(t, Mi) = (i/inun) x (a' - 1) x (M t / M^) 1 ^' (6) 

where w(t, Mi) is the weight assigned to a cell with age t 
and mass Mi, M max is the mass of the most massive cell 
in the simulation, f m ; n is the age of the youngest cell in 
the simulation and a 1 is the exponent of the mass func¬ 
tion. The weight is now a function of age and initial mass, 
where the most massive cell in the youngest bin has a 
weight of one, and corresponds to 1 cluster per log(f)- 
log(M) bin. This value is arbitrary and can change since 
later we will scale the total population to a different num¬ 
ber. When —a' = —2 is chosen, the weight depends on age 


and mass simply as: w(t, Mi) oc t/Mj. When the cells are 
then binned in age, mass or luminosity, the weights of each 
cell are counted and the bin values now represent number 
of clusters. The main advantage of using these weight as¬ 
signed cells spread equally in log(f) and log(M), is that it 
would otherwise be very time consuming to generate clus¬ 
ter populations with a well sampled mass function; the 
amount of clusters needed for that is too high. Also, the 
number of points per bin is now constant and it is straight¬ 
forward to create a variety of populations with different 
formation rates, di sruption time scales etc. in a short time. 
The model from lGieles et al.l 1 2005al) is extended in this 
work to create LFs. To this end the age and mass depen¬ 
dent magnitudes of clusters are taken from the GALEV 
SSP models based on a Salpeter IMF. When generating a 
LF, the total weight of the cluster population is scaled to 
a realistic number of clusters. Then only bins with weights 
greater or equal to one are kept. In this way the predicted 
maximum luminosity follows directly from the model, de¬ 
pending on the number of clusters (=total weight) and the 
input variables. 


5.1.2. Include stellar evolution and cluster disruption 


In a recent study ( Lamers et al 1l2005alf it was shown that 
there is a simple analytical description of the evolution 
of the mass of a single cluster as a function of time 
(see Eq. Q § 13.1.31) . It takes into account the effect of 
mass loss due to stellar evolution, based on the mass loss 
predicted by the GALEV models and it includes mass 
loss du e to tidal evap o ration based on the iV-body re¬ 
sults of lBaumgardt fc Makind l 2003ll . The mass as a func¬ 
tion of time according to these analytical formula agrees 
very well with the p redictions from IV-body simulations 
(lLamers et al . 2005alb We assume stars of all type are lost 
due to disruption. In reality more low mass stars will leave 
the cluster. This might slightly affect the magnitudes as a 
function of time. 


5.2. The luminosity function of a multi-age population 

In the first step to understand the LF we generate an an¬ 
alytic cluster population with log(f) ranging from 7 to 9. 
Cells are spread equally in log(f) and log(M) space and 
weights are assigned to correct for the logarithmic bin¬ 
ning in age and to create a mass function with —a 1 = 
—2.0. Cluster masses were assigned between 1O 3 M 0 and 
1O 12 M 0 . The upper mass was chosen unrealistic high, rep¬ 
resenting the lack of a physical upper limit to the most 
massive cluster in the mass function. The total weight 
is then scaled to 1000 clusters. In principle this step is 
not necessary, it tells us what the most luminous cluster 
in the sample is, since we will only use bins with values 
greater than one. (Note that due to the weights we use and 
the final scaling there are cells with weights much smaller 
than one.) Every cell is assigned a magnitude based on 
the GALEV models. The exponent of the resulting LF for 
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Table 1. Properties of the observed LFs. 


Galaxy 

—on 

— 02 

Break 

(mag) 

M* 

(mag) 

Foreground extinction (Ay) 
(mag) 

SMC 

-2.0 ±0.2 

- 

- 

- 

0.09 

LMC 

-1.9 ±0.1 

- 

- 

- 

0.13 

NGC 5236 

-2.1 ±0.1 

- 

- 

- 

0.218 

NGC 6946 

-1.7 ±0.5 

-2.4 ±0.1 

-8.9 ±0.4 

-10.2 ±0.6 

1.133 

M51 1 

-2.0 ±0.1 

-2.7 ±0.3 

-9.3 ±0.4 

-10.3 ±0.5 

- 

“Antennae” 2 

-1.9 ±0.1 

-2.8 ±0.4 

-10.3 ±0.4 

-11.4 ±0.4 

- 


1 Corrected for individual cluster extinction. 

2 Values from W99 for their sample of extended sources. 


this multi-age population is always —a = —2.0. Although 
clusters with log(f) = 9 have faded considerably, the final 
LF is still an addition of power-laws with exponents of 
— 2 , with no limit to the maximum mass and hence the 
maximum luminosity. For this reason, changes in the for¬ 
mation rate as well as bursts will have no effect on the 
slope. Mass dependent disruption will slightly change the 

slope fls isxn . 

5.3. The luminosity function without a maximum mass 

When the CIMF is not truncated, the mass of the most 
massive cluster observed will increase with increasing 
log(age/yr) due to the size-of-sample effect (§01 H03). The 
most massive cluster will depend on age as log(M max ) oc 
log(t), when —a' = —2. The dotted line in Fig. II 21 takes 
this increase of the mass and the fading of the most mas¬ 
sive cluster as a function of log(age/yr) into account. In 
two dex in age the maximum mass will also increase with 
two dex if the slope of the mass function is —2. This means 
that the most massive cluster at 10 9 years will be a factor 
of 100 more massive and hence 5 mag brighter than the 
one at 10 7 years. Since clusters fade about 4 mag between 
10 7 and 10 9 years, these two effects almost cancel out. Not 
only the mass of the most massive cell will increase, but 
also the second most massive etc. Therefore, the LFs at 
different ages will all have clusters sampled between the 
detection limit and more or less the same maximum lumi¬ 
nosity. So, the slope of the integrated luminosity function 
will be close to the slope of the CIMF for almost all lumi¬ 
nosities, which we found already in § 15.21 This effect can 
be seen in Fig. El for the LMC and the SMC clusters. At 
all ages we find clusters between the detection limit and 
more or less the same maximum luminosity. Integration 
over all ages will yield a LF which has about the same 
slope as the mass function. This has always been assumed 
to be logic, but here we show that this is only the case 
when there is no maximum mass. 

5.4. Truncation of the mass function 

L02 has already shown that when generating a multi-age 
population with a truncation at the high mass side, the LF 
gets steeper. This can be understood as follows: when the 
underlying mass function is truncated at some upper mass, 


the brightest cluster will be the most massive cluster in 
the youngest age bin. This holds when the MF is sampled 
up to the truncation limit in the youngest age bin. The 
older clusters have faded due to stellar evolution. Hence 
at the brightest end of the LF only the youngest clusters 
will contribute. Going to fainter magnitudes, more and 
more ages will contribute to the LF. At the luminosity of 
the most massive cluster in the oldest age bin and fainter, 
all ages will contribute, and therefore the slope faint-ward 
of that brightness will be the slope of the CIMF, similar 
to what we found in § 15.21 This is illustrated in Fig. El 
The decreasing line in the right hand panel represents the 
luminosity of a 1O 6 M 0 cluster between 10 7 and 10 9 years. 
The dark shaded region is where all clusters of all ages 
contribute to the same part of the LF. The light shaded 
region represent the clusters contributing to the bright end 
of the LF. 

When there is a physical mechanism that does not al¬ 
low clusters more massive than M max to be formed, there 
should be an observable bend at M\ = t max ). 

With our model we now generate a cluster population with 
ages between 10 7 and 3 x 10 9 yr and masses between 10 2 
and 10 6 and apply a detection limit at My = —7. The 
age range is based on the observed age range in M51. We 
excluded the clusters younger than log(age/yr)=7, since a 
large fract ion of these clusters will not survive the fir st 10 
Myr (e.g. lLada fe Ladal Efiol iBastian et al.ll2005bl). We 
apply no disruption or extinction. The total weight of this 
simulation is scaled to 1000 clusters. The brightest clus¬ 
ter now is the first bin at the bright side with a weight 
larger than 1 . The resulting LF is shown in Fig. El The 
faint part has an exponent of —a = — 2 . 0 , which is ex¬ 
pected, since the underlying mass function has an expo¬ 
nent of —o! = —2. On the bright side the slope is steeper 
(—2.5) and the break between the two power-laws is at 
My = —9.3 mag, which is close to the magnitude of a 
1O 6 M 0 cluster with an age of 3 x 10 9 yr ( My{t = 3 x 10 9 
yr; M = 1O 6 M 0 ) = —9.1 mag) . For NGC 6946 and M51 
the observed break occurs at My = —8.9 ± 0.4 mag and 
My = — 9.3 ± 0.4 mag respectively implies a maximum 
cluster mass of 6.9 x 1O 5 M 0 and 1.0 x 1O 6 M 0 . For com¬ 
parison, the bend in the LF of the “Antennae” galax¬ 
ies at My = —10.3 corresponds to a maximum mass of 
~ 2.5 x 1O 6 M 0 , assuming the oldest cluster in the sample 
is also 3 Gyr. We also fit a Schechter function (Eq.|5J) to the 
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Fig. 12. Right: Schematic picture of the two regions 
(dark and light grey) contributing to the LF when the 
underlying CIMF is truncated. The dark shaded region 
contributes to the faint end of the LF and the light shaded 
region to the bright end. The line indicates the brightness 
of a 10 6 Mg cluster. The dotted line indicates the bright¬ 
ness of the most massive cluster taking size of sampling 
effects into account. Left: Resulting LF. 


model as was done for the clusters in NGC 6946, M51 and 
the “Antennae” galaxies. When we fit this function to our 
model with a, IV* and M* as variables we get a very good 
fit (see bottom panel of Fig. 1 1311 . As was the case with 
the Schechter fit to the observed LF of NGC 6946 and 
M51, the bend now occurs brighter than when breaking 
the LF into two power-law distributions. When compar¬ 
ing the bend points from the Schechter functions to the 
model, we find maximum cluster masses of 4.4 x 10 5 Mq, 
4.8 x 10 5 M 0 and 1.3 x 10 6 M e for NGC 6946, M51 and 
the “Antennae” respectively. We note that the errors in 
the fitted break points and Schechter bends are large (~ 
0.5 mag), which corresponds to a factor 1.5 in mass. It is 
therefore hard to place a hard limit on the upper mass of 
the CIMF. 

5.5. Impact of different initial conditions 

Our assumed initial conditions of §EH might affect our 
results and in particularly the location of a break in the 
LF. We therefore created several cluster populations, with 
different initial conditions. 

5.5.1. Varying cluster formation rate 

For simplicity we have taken the cluster formation rate 
(CFR) to be constant over 3 Gyr in § 15.41 Especially 
when comparing with M51 and the “Antennae” galaxies 
this might be an over simplification. We generated mod¬ 
els with burs ts in the CFR, as is probably more realis¬ 
tic for M51 llGieles et al.l l2005a f. Bursts with increased 
CFRs of factors of 2-3 do not affect the location of the 
bend. For comparison with the LF of the “Antennae”, 



Fig. 13. LF of a synthetic cluster population with a CIMF 
with exponent -2, a mass function truncation at M max = 
10 6 Mq and a detection limit at My = —7. The top panel 
shows a double power-law fit, where the slope bright-ward 
of the bend is steeper than the underlying mass function. 
The bottom panel shows a fit with a Schechter function. 


we incr eased the CF R the last 160 Myr, based on the re¬ 
sults of I Wilson et al.l 1 2003lf . For a CFR which was 3/5/10 
times higher the last 160 Myr as before, the bend occurs 
0.3/1.5/1.8 mag brighter. This would make our observed 
maximum masses 1.3/4/5 times lower. As was shown in 
§ ?, the assumption of a constant CFR will always make 
it hard to directly compare the model with the data. 


5.5.2. The stellar IMF 

Our SSP models assume a Salpeter IMF. When a Kroupa 
was assumed, the corresponding mass derived from a bend 
in the modeled LF will be about a factor of 2 lower. To 
explain the bend location in M51, NGC 6946 and the 
“Antennae” galaxies, our derived maximum masses will 
also be a factor of 2 lower. 


5.5.3. The age range 

We assumed clusters in the age range of 10 7 < age 
< 3 x 10 9 years, based on the observed cluster ages in 
M51. When we would take 10 9 years as the maximum age 
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range, the location of the bend would be 0.8 mag brighter. 
The maximum masses of our observations would be a fac¬ 
tor of 2 lower in that case. If we would change the youngest 
cluster in our model from 10 7 years to 4 x 10 6 years, the 
modeled location of the bend does not change. However, 
when we assume that only 20 % of the youngest clusters 
(e.g . age < 10 7 years') survive t he gas removal phase (e.g. 
Fall 2004t IBastian et a,l.N20f)5hl) . it would place the bend 
0.4 mag fainter. That would make our maximum masses 
derived from observations a factor of 1.5 more massive. In 
addition, clus ters will probably loose 50-80% of their stars 
in this phase dKrouna Hr. HoilvIboO^l . which will make all 
our derived initial maximum masses a factor 2-5 more 
massive. 


5.5.4. The effect of disruption 

Applying cluster disruption according to Eq.Qaffects only 
the LF at the faint end. The low mass clusters are dis¬ 
rupted faster and that will result in a shallower mass and 
hence LF at higher ages. Since the young clusters also con¬ 
tribute at these brightnesses, it is hard to quantify this 
affect. When applying the short cluster disruption time of 
M51 (£4 = 10 8 years), we find an exponent of the LF on the 
faint side of a = 1.7 and the location of the bend does not 
change. It could explain why in the “Antennae” galaxies 
slopes significantly shallower than -2 are reported (W99) 
after various steps to removes stellar contamination. Also 
this value is close to the observed faint slope of the LF of 
NGC 6946 and it could explain why the slopes of the LF 
of the SMC and LMC are shallower than -2. 


5.6. Testing the bend scenario in different filters 


If the location of the bend in the LF is related to the most 
massive cluster in the oldest age bin (§ 15.411 . it should ap¬ 
pear at brighter luminosities in the red filters because old 
clusters have red colors. With this knowledge, we can test 
our theory on the cluster sample of M51, for which we have 
multiple filters available. With our model we can now pre¬ 
dict the LF and a possible bend in all filters. We used a 
constant formation rate and a disru ption time o f cluster s 
of £4 = 10 8 years, as was found bv ldeles et alJ ( 20f)5a b 
The CIMF was truncated at the high mass side. The clus¬ 
ters were generated and evolved according to method de¬ 
scribed in §EU The bend will occur at different magni¬ 
tudes for the different filters because the magnitude of the 
maximum cluster will be different in the different filters. 
The result is shown in Fig. El The top and bottom panel 


show the data in the F336W, F439W, F555W and F675W 
filters, roughly equal to U, B,V and R. In the top panel 
we plot the data and over plotted double power-law fits 
to the model. The break point between the two power- 
laws was based on the magnitudes of M max in these filters 
at 10 9 years according to the GALEV SSP models. The 
predicted bend is indicated in each band in the top panel. 


The break moves to brighter luminosities for redder 
filters as predicted, this is because clusters of 3 x 10 9 year 
are in general redder than younger clusters. For example, 
the B — V color at 10 9 year is 0.7 mag, so the bend in the 
E-band occurs 0.7 mag brighter than in the U-band. In 
the observed LF also the bend seems to move to brighter 
luminosities when going to redder filters. In the /7-band, 
the bend is predicted 0.8 mag brighter than the detection 
limit. We thus only see the steep part of the LF. The ob¬ 
served LF of the U -band seems also to be a single power- 
law, although a slight bend is visible just bright-ward of 
the detection limit. The steep part is much shallower in 
the F336W band than in the other filters due to the rapid 
fading in that filter. The data and the model agree very 
well on this. The bottom panel of Fig. 1141 shows the same 
data, but now the models were fitted with Schechter func¬ 
tions (Eq. 0. The fit to the model ends at the bright side 
at the last bin which had a value of 1 or higher. This means 
that this is the luminosity of the brightest cluster in the 
simulation. The brightest cluster in the model agrees well 
with the brightest cluster in the data. 


The bend in t he “Antennae” occu rs ~ 5 mag brighter 
in the Ks-band I Men gel et al ] 2 Qq 3) as in the E-banc 
(W99). The V — K model color of a cluster of 3 Gyr is 
around 3.2 mag. The E-band observations of W99 are not 
corrected for extinction, making the bend in E fain ter and 
most probably the observations by iMengel et al.l (2005) 
are hampered by blends, making the bend brighter in Ks. 


5.7. Effect of a bend in the mass function 

The observed bend in the LF of the “Antennae” galaxies 
ha s been i nter p reted a s a bend in the mass function (W99; 
lFritze-v. Alvenslebenlli 1999 1. With our model we can pre¬ 
dict how a bend in the CIMF would affect the LF. We 
generate a cluster sample between 10 ' and 10 9 years with 
a constant formation rate. We put a bend in the CIMF at 
M = 10 4 Mq, where the exponent is —a' = — 1 below the 
bend and —cd = —2 on the high mass end. The resulting 
LF is shown in Fig. El On the very faint side the LF has 
the same slope as the underlying MF below 1O 4 M 0 . Note 
that a slope of the LF of —a = — 1 shows up as a flat 
part in the magnitude distribution due to the conversion 
between the slope of the magnitude distribution and the 
slope of the LF (Eq. 0. On the bright side the LF has 
the same slope as the MF above 1O 4 M 0 (—a = —2). The 
transition region is caused by the fading of a 1O 4 M 0 clus¬ 
ter between 10 7 and 10 9 years. When we would fit two 
power-laws over the whole LF, the break would be some¬ 
where at My = —6 mag, i.e. much to faint to observe a 
bend in the MF at 1O 4 M 0 . 

5.8. Bend vs. no bend 

So if the bend in the LF (§ Ed can be explained with 
a truncated mass function, why do the other galaxies in 
our sample (SMC, LMC, and NGC 5236, § 14.1.114.l~^l 
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Fig. 14. The LF of 705 clusters with ages between 10 7 and 10 9 year in M51 in the F336W, F439W, F555W and F675W 
bands. Over plotted with solid lines are fits to modeled luminosity function, based on a clusters IMF with exponent 
-2.1, equal number of clusters and a truncation of the mass function at M max = 10 5 9 . Top: Fits to the model are two 
power-laws, with the break at M\ = t m ax )- The bend occurs to brighter magnitudes going to redder filters, 

since a cluster of age 10 9 years is brighter in the redder filters. Bottom: Same data, but now the fits to the model are 
Schechter functions. 



Fig. 15. Simulated LF of clusters which are formed with 
a constant formation rate between 10' and 10 9 years and 
a bend in the CIMF. Below 10 4 Mq the mass function has 
an exponent of —a! = — 1 and above 10 4 M© the exponent 
is —a' = —2. The mass function has no truncation at high 
masses. At the bright side the LF has the a slope of a = 2 
and on the very faint side the slope is a = 1. 

not show this bend? The answer is sampling statistics. If 
the bend is a physical feature occurring at a given mass, 
then it should always occur at a predictable luminosity 
(assuming the age distribution is known). Of the galaxies 
studied here, NGC 6946 and M51 (and the “Antennae”) 
are the only ones where the LF is sampled above My = 


— 10. In the SMC, LMC and NGC 5236 we do not sample 
the LF up to such bright magnitudes, so if the truncation 
occurs at the same mass as in N6946/M51 we would not 
be able to see it. 

6. The maximum cluster luminosity 

Now that we have shown that there are arguments to 
believe that the CIMF is truncated, it is interesting to 
compare how this might a ffect t he m aximum cluster lu¬ 
minosity in a galaxy. Kviiitmoril li2003h plotted the maxi¬ 
mum cluster luminosity in a set of galaxies is a function 
of the number of clusters above My = —9 mag. A fit to 
the data is close to the expected relation based on the 
size-of-sample effect and an exponent of the luminosity 
function of a = 2. L02 has proposed a relation between 
the maximum cluster luminosity and the product of the 
observed area (A) and the area-normalized SFR (£sfr) 
of the galaxy, assuming that the cluster formation rate is 
proportional to the SFR. The increase of the maximum 
cluster luminosity with the product A*SFR is consistent 
with the observed increase. 

In the top panel of Fig. nnwe show the brightest cluster 
as a function of log (N(M y < —8.5)) from the sample of 
lLar gen Richtlerl feoQflf with filled circles_ajidj^hitm 2 rd 
ll2003lf with open circles. The sample of lWhitmord 1120031 
has been shifted to the right with 0.2 dex, which is the 
expected number of clusters missed between My = —8.5 
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and My = —9.0 based on a luminosity function with expo¬ 
nent a = 2. The A log(TV) in a certain magnitude interval 
relates to the slope of the LF as 


Alog[7V(AM)] = 0.4 (a-1) AM, 


(7) 


where AM is the magnitude interval and a is the slope of 
the LF. We h ave updated t he data for NGC 1569, based 
on the data oflH unter et~ahl I 2000h . This was an outlying 
point injwhitmore ( 20031 1. but we found more clusters in 
iHunter et ali I200d). 


We have over-plotted the expected increase in max¬ 
imum cluster luminosity with log(N) assuming the lu¬ 
minosity function is a single power-law distribution and 
based on a range of exponent (1.7 < a < 2.7), which 
is the observed range of slopes in this work. The fit to 
the closed circles is shown with the solid line, which falls 
nicely within the predicted area and corresponds to a ex¬ 
ponent of the LF of a = 2.45 ± 0.30. Note that although 
this value is slightly higher than found in § ro it is still 
within the ler error. Also we here use an indirect method 
to determine the slope of the LF, based on a small number 
of data points. 

L02 has quantified the expected scatter in the max¬ 
imum cluster luminosity. The ler deviation is 1.04 mag 
from the mean and is independent of the number. We 
over-plotted the 3er deviations. To show the effect of the 
expected scatter, we show a simulated set of maximum 
cluster luminosities as a function of \og(N) in the middle 
and bottom panels of Fig. ED The input exponent of the 
simulated LF was the value found from the fit to the obser¬ 
vations (2.45). The same number of simulated maximum 
cluster luminosities as was observed is shown in the mid¬ 
dle panel and a larger number in the bottom panel. The 
simulated luminosities show that luminosities close to the 
3cr error line are expected even with the low number of 
galaxies. The 3cr scatter is a little bit higher on the bright 
side of the fit, which is because of the power-law nature 
of the LF function. 

For the galaxies studied here, we find that the lumi¬ 
nosity of the brightest cluster can be accounted for by 
sampling statistics. This is in contrast with the most mas¬ 
sive cluster. This can be understood by looking at Fig. |2] 
and Fig. El There we can see that the most luminous clus¬ 
ter in each galaxy is always a young cluster (log(age/yr) 
~ 7). The most massive cluster, however, is usually in the 
older bins, since these cover a larger time in which it is 
more likely to form a more massive clusters. The physical 
upper limits we have derived in this work are not reached 
when sampling a CIMF during 10 7 years. 


7. Discussion 

A truncation of the cluster mass function raises the ques¬ 
tion what physical mechanism could be responsible for this 
cut-off. An answer might be present in the mass spectru m 
of giant molecular clouds (GMCs). IWilliams fc McKeel 
( 1997 1 have shown that there is a distinct cut-off at the 



Fig. 16. The maximum cluster luminosity in a galaxy vs. 
the number of clusters above a certain magnitude limit. 
Top: Observed maximum l uminositie s for 37 galaxies, 
compiled f rom a sample by Whit morel ll 20031) with open 
circles and lLarsen fc RichtleJ ifeoQOfl shown in full circles. 
The open circles are shifted 0.2 dex, since that sample 
is for clusters with My < —9 mag, while the sample of 
lLarsen fc Richtleil ( 2000 1 is for clusters with My < —8.5. 
The grey area represents the expected range based on 
the size-of-sample effect and the observed range of a in 
this study (—2.8 < —a < —1.8). The full line is a fit to 
the data. Some special cases are indicated as: M51 (o), 
NGC 1569 (A). The middle and bottom panel show 
simulated maximum cluster luminosities for similar num¬ 
ber of galaxies as the observations and a large number 
respectively. The simulations were based on —a = —2.45. 
In all three plots the 3<r deviation predicted by L02 is 
over-plotted with a dashed line. 


high mass end of the mass function of GMCs in the Milky 
Way. There are no clouds observed above 6 x 1O 6 M 0 , while 
there are ~100 GMCs expected above that cut-off, based 
on statistical arguments. This might be caused by shear 
effects which do not allow clouds to grow larger than a 
certain size. The size and the mass of clouds are directly 
related by the virial theorem. Assuming that star clusters 
are formed with a star formation efficiency of e = 0 . 01 — 0.1 
and that one single GMC produces one cluster, this im¬ 
plies it is not possible at this moment to form clusters 
more massive than 0.6 — 6.0 x 1 O 5 M 0 in the Milky Way. 
The Milky Way has formed a few massive clusters in the 
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last few My r, like Westerlund 1 with a mass of ~ 1O 5 M 0 
IClark et al. '200 5 • and the Arches and Quintuplet clus- 
ters i n the Galact ic center with masses of around 1O 4 M 0 
llFiger et al.lll999 :i ) . 

A rough estimate can show that we would have ex¬ 
pected more massive objects. Westerlund 1 is the most 
massive young cluster in the Milky Way disk. Of course 
we have to look at the size-of-sample that was used to dis- 
cover it. Wes terlund 1 is about 4.5 kpc away from the sun 
( Clark et al i ! 2005f ). Based on the most rece nt cat alog of 
open cl usters by iKharchenko et al | (120051) . lLamers et al.l 
~ 2005all estimated the masses of clusters in the solar neigh¬ 


borhood. Within 600 pc from the sun, the catalog is com¬ 
plete until 200 Mg, where about 100 clusters are found 
above that mass. The initial mass of the most massive clus¬ 
ter observed is estimated to be 3 x 1O 4 M 0 . The size of sam¬ 
ple predicts M max = N M min = 100*200M e = 2xlO 4 M 0 . 
So within the solar neighborhood the mass of most mas¬ 
sive cluster is determined by the size-of-sample effect. 
Assuming that the number density of clusters is constant 
until Westerlund 1, we can predict the expected maximum 
mass, based on statistical grounds. Based on the distance 
of 4.5 kpc, the expected number of clusters above 200 M 0 
is about (4500 2 /600 2 = 56) times larger. This implies we 
would expect a maximum mass of M max = 1.5 x 1O 6 M 0 , 
Although Westerlund 1 is a very massive cluster with 

10 5 M 0 , we would have expected at least one cluster above 

10 6 M & . Given that Westerlund 1 is highly obscured, it is 
unlikely that a 10 times more massive object would not 
have been found yet within a distance on 4.5 kpc from 
the sun. In fact, the mass of Westerlund 1 might be a rea¬ 
sonable uppe r limit, given that the mo st massive GMC is 
6 x 1O 6 M 0 ( Williams fc McKee! i'l997 f and a reasonable 
star formation efficiency of a few percent. 

It is not clear what causes the cut-off of the cloud 
mass function. The maximum cluster mass seems to be 
environment dependent, since in the “Antennae” galax¬ 
ies clusters of a few times 1O 6 M 0 are being formed. 
This could be because the maximum GMC mass is 

or because the star forma- 


Tils 


aU 2C 


higher 

tion efficiency is much higher. The star formation ef¬ 
ficiency is e xpected to be hi gher in high pressure en¬ 
vironments llElmegreen fc Efremov! Il997l) . which is the 
case in shock induced star formation environments like 
galaxy mergers 


). The super massive clus¬ 
ter W3 in NGC 7252, with a dynamical mass of 1O S M 0 
llMaraston et al.l)2004fl . exceeds the limits we observe by 
almost two orders of magnitude. It was probably formed in 
the merger process of two gas-rich spirals and seems to be 
the ti p of a continues power-law distribution I Miller et ID 

1M7). 

As a comparison, the “Antennae” galaxies have not 
merged their nuc lei vet. S imulations of the m erging of two 
gas-rich spirals ( Barnes fc Hernouistl ll991 1 have shown 
that in the third encounter the gas disks merge with rel¬ 
ative velocities of more than 500 kins -1 . Perhaps more 
massive clusters can be formed in these merging stages, 
due to the merging of clusters. Clusters also seem to form 


in complexes of multiple clusters and stars, w ith simi¬ 
lar pr o perties of their progenitor GMCs i Ba.stian et ID 
). iFellhauer fc Krmrnal ( 2005 1 have shown by nu¬ 


merical simulations that these complexes can collapse and 
form a single, diffuse, ultra-nrassive object, which might 
be the way W3 and other very massive clusters can be 
formed. During the merger shear-free regions exist, like 
in the overlap-region in the “Antennae” galaxies, where 
GMCs can grow bigger and more massive than in the 
galaxies we have studied here. 


8. Conclusion 

We have compared observed star cluster luminosity func¬ 
tion in five galaxies with analytical cluster population 
models. Our main results can be summarized as follows: 

— there are no clusters in M51 more massive than M ~ 
1 x 1O 6 M 0 , although they are predicted by the size-of- 
sample effect. When comparing the maximum cluster 
mass in increasing log(age/yr) bins, the LMC and SMC 
cluster population show an increase consistent with the 
size-of-sample of effect. The cluster population of M51, 
however, shows a much shallower increase. This sug¬ 
gests a physical upper limit to the masses of clusters 
M51, although the shallow increase can also be repro¬ 
duced by a combined effect of cluster disruption, infant 
mortality and an increasing cluster formation rate. 

— when comparing the xZ results of different function 
fits to the five galaxies in our sample, we find that the 
LF of the SMC, the LMC and NGC 5236 can be well 
approximated by a power-law (N dL oc L~ a dL), with 
1.9 < a < 2.1, while the LF of NGC 6946 and M51 are 
slightly better approximated with a double power-law 
or Schechter function. 

— when fitting a double power-law function to the LF of 
NGC 6946 we find a break point at My = —8.9 ± 0.4 
mag. Faint-ward of the bend a power-law with expo¬ 
nent — 1.7 ± 0.5 can be fitted. Bright-ward of the bend 
an exponent of —2.4 ± 0.1 is found. The LF can also 
be well fitted by fit a Schechter function with a bend 
at My = —10.3 ± 0.5 mag. 

— the LF of M51 shows a break at My = —9.3 ± 0.4 
mag. Faint-ward of the bend a power-law with expo¬ 
nent —2.0 ± 0.1 can be fitted. Bright-ward of the bend 
an exponent of —2.7 ± 0.3 is found. The LF of M51 is 
also well fitted by a Schechter function with a bend at 
My = —10.3 ± 0.5 mag. 

— the cluster LFs can be reproduced with a synthetic 
cluster population model. The bend in the LF of 
NGC 6946, M51 and “Antennae” galaxies can be ex¬ 
plained with a truncation of the cluster mass function 
at M max = 0.5 - 1 x 1O 6 M 0 (M51/NGC 6946) and 
1.3 - 2.5 x 1O 6 M 0 (“Antennae”). 

In a follow-up study JCieles et al.l2006bil we present an 
improved LF of star clusters in M51 based on HST/ACS 
data, taken as part of the Hubble Heritage project. 
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